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We investigate the fermion determinant of a new action on a (1 + D)-dimensional lattice for 
SU(2) gauge groups. This action possesses the discrete chiral symmetry and provides 2 D -component 
fermions. We also comment on the numerical results on fermion determinants in the (1 + D)- 
dimensional SU(3) gauge fields. 
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I. INTRODUCTION 

o ■ 

As is well known the lattice formulation of fermions has extra physical particles or breaks the chiral symmetry. 
£NJ | This is unavoidable under a few plausible assumptions . Several methods have been proposed to deal with this 
difficulty. Wilson's formulation ||, which is one of the most popular schemes, eliminates the unwanted particles 
with an additional term which vanishes in the naive continuum limit. However, this formulation sacrifices the chiral 
symmetry. An alternative scheme is the staggered fermion formulation proposed by Kogut and Susskind This 
\q ' scheme preserves the discrete chiral symmetry and in this point the staggered fermion has an advantage over the Wilson 
fermion. But, the staggered formulation describes a theory with 2~ ~ degenerate quark flavours (2 1+D components) 
in (1 + D) dimensions, while there is no restriction on the flavour number in the Wilson formulation. Recently, it has 
been shown that lattice fermionic actions with the Ginsparg- Wilson relation Q have an exact chiral symmetry and 
are free from restriction on the flavour number. But, these actions cannot be "ultralocal" which makes numerical 
(*— ^ | simulations complicated. 

In the recent papers we proposed a new type of fermionic action on a (1 + Z3)-dimensional lattice. The action 
is ultralocal and has discrete chiral symmetry. On the Euclidean lattice the minimal number of fermion components 
is 2 D , which should be compared with 2 1+D of the staggered fermion. When dynamical fermions are included, the 
numerical feasibility relies on the reality and positivity of the fermion determinant. In the previous paper [g| we 
investigated, analytically and numerically, the fermion determinant of our new action in the (1 + l)-dimcnsional U(l) 
lattice gauge theory. We showed the reality of our fermion determinant under the condition fixing the global phase of 
link variables along the temporal direction. By a similar discussion to the U(l) gauge group, we could also find the 
reality and the positivity of our fermion determinant in the (1 + l)-dimensional SU(N) lattice gauge theory. 

In this paper we analytically show that our fermion determinant with the SU(2) gauge fields is real and positive in 
• (1 + D) dimensions. We also comment on the numerical results of the fermion determinant in the (1 + _D)-dimensional 
SU(3) gauge fields, and discuss the effectiveness of our new action for SU(2) and SU(3) lattice gauge theories. 

II. NEW FERMIONIC ACTION 

CZ \ 

In the recent paper [Q, we proposed a new fermionic action on the Euclidean lattice. Though the action keeps the 
discrete chiral symmetry like the staggered fermion action, the fermion field has 2 D components in (1 + D) dimensions. 
In this section we briefly sketch our formalism for later convenience. 

The action can be written with a fermion matrix A as 

Sf = ^V4 A nmV'm, (2-1) 

where the summation is over lattice points and spinor indices, and our fermion matrix is defined by 

A = 1 - SlU E . (2.2) 
Here Ue is the Euclidean time evolution operator and is the unit shift operator defined by 

S^(x°, x 1 , . . . , a?, . . . , x D ) = ifj(x°, x 1 , . . . , x" + 1, . . . , x D ) (/x = 0, 1, . . . , D). (2.3) 
We require that the propagator has no extra poles and find that Ue has the form 
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{ iX * - S l) + (! - Y >) (Si-2 + fit) } , (2.4) 

where is the ratio of the temporal lattice constant to the spatial one. The spinor matrices X's and y's should 
satisfy the following algebra: 

2 

{Xi,Xj} = <5y, 

{^,^} - 0, (2.5) 



where i and j run from 1 to £>. The matrix 2(5ij/rE + 1) is positive definite for any positive te, therefore X's and 
y's can be assumed hermitian, 

x\=X h Y}=Yi, (2.6) 
The matrices X's and Y's can be expressed by the Clifford algebra: 

r+=r n , {r„,r m } = 2S nm ( n ,m = i,...,2D) (2.7) 

in several ways. One is 

— D 

-Ti, Yi=J2 a io T D+j, (2.8) 



X, , 



as was used in the previous paper. Another one is 



X, = J — Ti, Y = J —T D+l + T 2D+1 , (2.9) 
\ r E V te 

where T2D+1 is 

T2D+1 = (-i) D TiT 2 ■ ■ - T 2 d- (2.10) 

The latter is more convenient than the former for later use. The dimension of the irreducible representation for T's 
is 2 D and accordingly ijj has 2 D components. 

The interaction of the fermion with gauge fields is introduced by replacing the unit shift operators by covariant 
ones: 

-> S M (x) = U XtX+(i S^, (2-11) 

where fi is the unit vector along the /i'th direction, and U x ^ y is a link variable connecting sites x and y. 
The fermion matrix Eq.(2.2) and the time evolution operator Eq.(2.4) become 

A(a;) = 1 - Sl(x)U E (x), (2.12) 

and 

D 

u E (x) = 1 - E y { lXl ( Sl{x) - s ^ x) ) + (1 " Yi) ( Sl{x) - 2 + s * {x) ) } • (2 - 13) 
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III. FERMION DETERMINANT FOR SU(2) CASE 



In this section we analytically study the determinant of our fermion matrix in SU(2) gauge fields. 
First, in the (1 + l)-dimensional case, the complex conjugation of U E (x) is 

U* E {x) = 1 - Y {~ iX * { S *W + ^'W) + (! - { S *( x ) ~ 2 + S *i( x )) } > 
where we can write 

3 

Sx(x) = ao(x)l +i s ^ j a. i {x)T i , (3.2) 
i=i 

since the link variables in Si(x) are SU(2) gauge group elements. Here ao(x) and cti(x) are real and depend on lattice 
points and T1.2.3 are the Pauli-matrices: 

ti = ( 1 ) ' T2 = ( °i 7 ) ' rs = ( -1 ) • (3 - 3) 

Then we have 

S 1 (x)t 2 =t 2 S* 1 (x). (3.4) 
By the same discussion for other unit shift operators Sq and S\ , we also have 

Sl(x)r 2 = r 2 S^(x), S\(x)t 2 = r 2 Sf(x). (3.5) 
If we can find the matrix T such that 



XfT = -TXi, 
Y{T = TYt, 



(3.6) 



it is easily shown that 

(r t 2 ) (sl{x)u E {x)) * (r ® t 2 ) = sl(x)u E (x). (3.7) 

For example, we make the following choice: 

X! = J —n, Y x - ./ — t 2 + r 3 , (3.8) 
V r E V r E 

the matrix T acting on two components fermi fields defined by 

T = r 3 (3.9) 



satisfies Eq.( 3^ ). The Eq.( p.7| ) implies that if A is some eigenvalue of our fermion matrix A(x), then A* is an eigenvalue 
of A*(x) and thus also of A(x). Therefore eigenvalues of A(x) are either real or come in complex conjugate pairs. 
From the above discussion we can prove the reality of our fermion determinant for the SU(2) gauge groups. 
Next we show its positivity. We define 

r' = (r®T 2 )x, (3.10) 

where K is complex-conjugation operator. We find 

T'r' - (r 3 ® r 2 ) K (r 3 ® r 2 ) K = (r 3 ® r 2 ) (r 3 ® (-r 2 )) = -1, (3.11) 



and from the relation Eq.(3.7) we can show 

[T',A(x)] = 0. (3.12) 
For a real eigenvalue A# of A(a;) and the eigenvector vr for this eigenvalue, from Eq.( 3.12] ) we obtain 
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A(x)T'v R = T'A(x)v R = T'XrVr = XrT'vr. (3.13) 
Suppose T'vr = cvr, then we find 

T' 2 v R = T'cvr = c*T'v R = \c\ 2 v R , (3.14) 



which is inconsistent with Eg. ( 3. 11 ), so that T'vr is different eigenvector for the same eigenvalue. Therefore the 
eigenvalues on real axis are degenerate in pairs and the determinant of A(x) is positive. 

The above proof of the positivity of our fermion determinant for the SU(2) group can be expanded to higher 
dimensions. We can make a fundamental representation for the Clifford algebra with 2D elements T n (n = 1, . . . , 2D) 
using direct products of the Pauli-matrices: 

Ti = T\ ® T3 <g> • • • ® T3 

D-l 

r 2 = l <g> n (g) r 3 <8> • • • <g> t 3 

D-2 



T D = 1 ® ■ ■ • ® 1 ®Tl 

D-l 

Td+i = t 2 ® r 3 ® • • • ® r 3 



(3.15) 




r>-2 



r 0+i = 1 ® • • • ® 1 ®r 2 ® r 3 <8> • • • <g> r 3 



^2D = 

where i runs from 1 to D. It can be easily checked that r n 's satisfy the relation 

{r„,r m } = 2<5 m „ (n,m = l,...,2Z>). (3.16) 
Moreover, we can see that the matrix Ti is real and the matrix T^+i is pure imaginary: 

rj = r<, T* D+i = -T D+i (i = l,...,D). (3.17) 
From the anti-commutation relation, we find the hermite matrix T2_d+i which anti-commutes with all r„'s: 

T2D+1 = (— i) D TiT2 ■ ■ - T 2 d 

= T3®T3<g>---®T3. (3.18) 

Clearly, the matrix T2_d+i is hcrmitian and the square of this matrix is equal to the unit matrix. Thus, we have 

r 2 D+ir*r 2 _D+i = -r i; ^iqa 

2D+1L D+i L 2D+1 — J-D+i- 



In (1 + D) dimensions, the relation Eq.(3.6) is rewritten as follows: 

x*v = -rx i: 
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FIG. 1. The spectrum of A(x) in the complex plane on a 4 x 4 x 4 lattice in SU(2) gauge group. 



Then, for the representation of Eq. (|2.9|) 



we find 



Xj, — \ l — Ti 



Yi — \ i — Tp+i + r2D+i, 



223+1- 



(3.21) 



(3.22) 



Since the eigenvalues of A(x) always consist of complex conjugate pairs and degenerated ones on real axis, we conclude 
the determinant of our fermion matrix is positive in SU(2) gauge fields in any dimensions. 

Now we show a numerical evidence. Fig.[l] shows the spectrum of our fermion matrix in a typical background 
configuration of link variables for SU(2) gauge group in (1 + 2) dimensions. We find that their distribution is 
symmetric with respect to the real axis as expected. Similarly in (1 + 3) dimensions we can numerically confirm the 
symmetry with respect to the real axis, and the positivity of the determinant. 



IV. DISCUSSION AND SUMMARY 



In the previous paper [gj we reported analytical and numerical results on the fermion determinant of our new action 
in (1 + 1) dimensions. In the case of U(l) gauge group, we were faced with the problem of convergence in numerical 
simulations. The cause of the poorness of the convergence is that the summation of the det(l — Sq(x)U e {x)) = 
det(l - e %e Sl(x)U E (x)) over arbitrary phase angle is canceled out accidentally. The element e comes from the 



one sided time difference operator Sq(x) with (x) replaced by 9q{x) + 0, i.e. Sq(x) — e lB Sl(x), where 9q(x) is 
defined by U x = e l6,(l(x ) ||. Therefore we must have control of the phase angle in order to get good convergence 
in the (1 + l)-dimensional U(l) gauge theory. In fact we analytically showed that our fermion determinant is real for 
all configurations and positive for most configurations under the T-condition (9q(x) = 0), which corresponds to the 
temporal gauge condition on the infinite lattice, or the GT-condition (y~] r 0q(x) = nit : n = even), which is achieved 
by a gauge transformation on the infinite lattice. It was also verified numerically. On the other hand we got good 
convergence without any conditions in (1 + l)-dimensional SU(N) case, because the element like e does not belong 
to the SU(N) group. 
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(a) SU(2) 



(b) SU(3) 
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FIG. 2. The distribution in the complex plane of our fermion determinants in a) SU(2) and b) SU(3) for each configurations 
of 2000 Monte Carlo iterations after getting good equilibrium, i.e. after 2000 iterations, on a 4 x 4 x 4 lattice at j3 = 2.0. 



The above discussion is applicable to higher dimensions to a certain extent. In SU(2) group our fermion determinant 
is analytically shown real and positive in any dimensions. In Fig.^|(a) we give the numerical evidence that the 
determinant is real and positive in (1 + 2) dimensions. In the case of SU(3) group, we cannot prove the reality of 
the determinant. But from Fig.§(b) we see that the distribution of the determinant is concentrated near the real axis 
without any conditions and the phase angle of the determinant is small. We have obtained similar results in (1 + 3) 
dimensions. When the phase angle of the fermion determinant is small enough, we can neglect the phase factor and 
make use of |detA(a;) instead of detA(x). In the above numerical simulations link variables are updated by the 
Metropolis method and determinants are calculated by the LU decomposition. So there are no systematic errors in 
the determinants. 

In conclusion, we believe that our new fermionic action is a profitable formulation for the numerical simulations of 
SU(2) and SU(3) lattice gauge theory. 
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